An algorithm for characterizing SLOCC classes of multiparticle entanglement 
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It is well known that the classification of pure multiparticle entangled states according to stochastic 
local operations leads to a natural classification of mixed states in terms of convex sets. We present 
a simple algorithmic procedure to prove that a quantum state lies within a given convex set. Our 
algorithm generalizes a recent algorithm for proving separability of quantum states [J. Barreiro et 
al., Nature Phys. 6, 943 (2010)]. We give several examples which show the wide applicability of our 
approach. We also propose a procedure to determine a vicinity of a given quantum state which still 
belongs to the considered convex set. 

PACS numbers: 03.67.-a, 03.67.Mn, 03.65.Aa 



I. INTRODUCTION 

The importance of quantum entanglement for quan- 
tum computation as well as for many other applications 
in quantum information processing has raised many fun- 
damental questions regarding its characterization [l| . In 
mathematical terms, a quantum state is said to be sep- 
arable, if it can be written as a mixture of projectors 
onto product states, otherwise it is entangled. Much 
work has been devoted to the development of criteria 
(in particular, the well-known tool of entanglement wit- 
nesses) which can prove that a quantum state is entan- 
gled, which means that it is outside of the convex set of 
separable states [2[. Interestingly, methods which prove 
that a mixed quantum state is within the set of separable 
states (e.g., by providing an explicit decomposition into 
product states), are less well known. Nevertheless, for 
some cases explicit decompositions are known [3|-|6[ and 
recently even some algorithms for this task have been 
developed 0-[il- 

If more than two particles are considered, the problem 
becomes more complicated, since different classes of mul- 
tiparticle entanglement exist. One possibility uses the 
notion of stochastic local operations and classical com- 
munication (SLOCC) [iflllll- For this notion, one can 
again ask, whether a given state can be decomposed into 
states of the same SLOCC class, which leads for the case 
of three spin-1/2 particles to the well known classification 
into GHZ- and W-states [I4I . Distinguishing these classes 
is a hard task, for partial results see Refs. [12, Il3- Es- 
pecially if one wishes to prove that a given state is within 
an entanglement class (such as the W-class), no general 
methods are known. 

In this paper we propose an algorithm which allows 



'Electronic address: kanipermann@fliphy. uni-duesseldorf.de] 

TElectronic address: otfried.guehne@ uni-siegen.de| 

-f Electronic address: wilmott@fi.mu ni.czl 

§ Electronic address: |brussQthphy7uiii-duesseldorf.de| 



to prove that a given mixed state belongs to a specific 
SLOCC entanglement class i.e., a decomposition exists 
where the pure states belong to the specified entangle- 
ment class. During the iterative procedure pure states 
and probabilities of the decomposition are determined as 
well as the "rest" operator which has an increasing mixed- 
ness during the iterations. In case of convergence we 
determine an explicit decomposition of the initial mixed 
state in terms of a convex combination of projectors onto 
pure states with the desired properties as well as a "rest" 
which is verified to be fully separable. It is not known 
whether the algorithm converges in all cases, but the 
method is easy to implement and it turns out to be well 
working in practice. Our algorithm is a generalization of 
the algorithm for proving separability from Ref . [9[ . 

Additionally we present a simple method to determine 
a lower bound on the e-ball of the calculated decomposi- 
tion [ij. This we achieve by constructing a cross poly- 
tope inside the convex set spanned by our decomposition. 
We find a ball with respect to the Hilbert-Schmidt dis- 
tance for states in this convex space. This method helps 
to verify that such properties also hold in case of reason- 
able small "experimental" errors. In some cases it also 
helps to get an idea of how far the state is at least away 
from sets with different properties. 

The paper is organized as follows. After an overview 
about multipartite entanglement and SLOCC entangle- 
ment classes in Section |ll] we present in Section Hill an it- 
erative algorithmic procedure to incrementally determine 
constituents of a decomposition of a given mixed quan- 
tum state. The main ingredient, the maximization of the 
overlap of a pure state with a given density operator un- 
der SLOCC operations, will be presented in Section IIVI 
and the properties of the algorithm in terms of conver- 
gence as well as scaling is discussed in Sections FVllVHI 
It follows in Section IVIIII the procedure of determining 
an £-ball. At the end a scries of examples, where de- 
compositions are determined and e-balls are calculated, 
is presented in Section IIXI and compared with results al- 
ready known in the literature. We conclude with possible 
further improvements and limitations of our procedure. 



II. DEFINITIONS AND NOTATIONS 

Let us first consider pure states. Generally, a pure 
composite quantum state of n particles is called entangled 
if it cannot be written as a tensor product of local states 
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where the |^i) are states on a subset of all n particles. De- 
pending on the value of k one can further distinguish the 
number of parties involved which arc not of tensor prod- 
uct form. Consequently, the states are said to be bisepa- 
rable (A; = 2), triseparable (fc = 3), up to n-separablc. If 
a state does not possess any tensor product structure it 
is called genuine multipartite entangled. 

This classification of pure state entanglement can be 
refined to the equivalence under stochastic local oper- 
ations and classical communication (SLOCC) [l3, KLlj . 
Physically speaking, SLOCC operations can be imple- 
mented with nonzero probability via local operations and 
classical communication, i.e. a single copy of |'0) can be 
mapped onto \4>) using local operations with probability 
p > 0, but with probability (1 —p) some other state may 
result. If two states can be converted into each other 
via SLOCC, this implies that both states are in princi- 
ple useful for the same tasks in information processing, 
albeit the efficiency might be different. We denote such 
a class of SLOCC equivalent states by C. 

Mathematically speaking, a general SLOCC operation 
can be represented by the action of local operators i.e., 
^SLOCC = (8); ^ij where Ai are arbitrary operators act- 
ing on the ith party. An SLOCC operation maps the 
initial state \ip) to \(j)) by 
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where Af denotes the normalization. 

For mixed quantum states shared between n parties a 
state is called entangled if it cannot be written as a con- 
vex combination of an n-fold tensor product of projectors 
onto pure states [l^ i.e.. 
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One can extend this naturally by considering fc-separable 
states. Finally, a state is genuine multipartite entangled, 
if it cannot be written as a mixture of biseparable states. 
For three qubits and pure genuine entangled states 
there exist two types of entanglement classes which are 
not SLOCC equivalent [l3|: The two representatives are 
the GHZ state and the W state 
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(|000) + |111)), 

(|001) + |010) + |100)). 



(4) 



Any pure entangled state can either be transformed into 
\GHZ) or \W), but these two states cannot be converted 
into each other. The analysis and the hierarchy of the set 
of mixed W/GHZ states was then developed in Ref. jl3 |. 
For more than three qubits already an infinite number of 
inequivalent SLOCC classes exist |lll |. 



III. THE MAIN IDEA FOR THE ALGORITHM 

A. Structure of the problem 

In this section, we will describe the main idea from 
Ref. [9| to design an algorithm for proving that a quan- 
tum state g belongs to a given class C. In the most gen- 
eral case, the task is to prove that a quantum state p is a 
convex combination of some projectors onto pure states 
\(t>k)- This means that we can write 



Q^ ^Pk\<l>k){<l>k 



(5) 



where the pk form a probability distribution. The |(/>fc) G 
C are states within a specific class C: For instance, if 
one wishes to prove that g is fully separable, then C is 
the class of pure fully product states, or if g should be 
proven to belong to the W-class for three qubits, then C 
is the class of pure W class states (that is, the SLOCC 
orbit of \W)). In the following, we will denote the set of 
density matrices that can be decomposed as in Eq. (O as 
conv(C), meaning the convex hull of C 

In order to design an algorithm to check whether g can 
be decomposed as in Eq. ([5]) we will use the following two 
facts: 

(a) Convexity: First, the set of density matrices with 
a decomposition as in Eq. ([5|) forms a convex set i.e., if 
ga and gt are in conv(C), then (1 — p)ga + pgt is also 
in conv(C). This is indeed obvious by definition and it 
will be used as follows: Assume that we have three states 
ga, gt, and gc which obey 

Qb = - (Qa-PQc) ^ Qa = {^-p)Qb+PQc (6) 

1-p 

where < p < 1 and gc lies in the set conv(C). In 
this situation, if we can prove that gi, S conv(C), then 
ga G conv(C) must hold, too. We will use this fact in 
terms of the first part of Eq. ([6]): starting from ga we 
will subtract a tciuipgc with gc = \4)k){4>k\ and \(j}k) £ C. 
If we can then show that g^ G conv(C); this implies that 
ga e conv(C). 

(b) Highly mixed states are in conv(C).' As a second 
fact we need statements which imply that highly mixed 
states are in conv(C). This, of course, requires a speci- 
fication of the degree of mixedness and depends on the 
structure of C. For instance, if we consider a bipartite 
N X M system and if C denotes the set of pure product 
states, then it has been shown that if 
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then g is separable, that is g € conv(C) [1^1 ■ Similar 
results have been obtained for other situations |17H20l | 
For instance, an A^-qubit state with A^ > 3 for which 



tr(gn < 
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holds, is fully separable [2l|. This condition will be used 
as a termination condition in the algorithm below. 

It should be noted that there are cases, where a con- 
dition as Eq. (I?!) is not directly given. For instance, if C 
is the set of all symmetric product states, then the max- 
imally mixed state is clearly not in conv(C), as conv(C) 
consists of matrices acting on the symmetric space only. 
Even for the identity operator I5 on the symmetric 
space, a condition like Eq. ([7]) is not straightforward to 
derive, since there are non-symmetric states close to Is 
which are not in conv(C). In our paper, such problems do 
not play a role, a detailed discussion of symmetric states 
will be given elsewhere. 



B. The algorithm 

Now we can formulate the iterative algorithm to prove 
that a state g is in conv(C). The algorithm consists of 
the following steps: 

1. Take the input state g as gk with A: = 1. 

2. Consider the optimization problem 



max 

I0)ec 



\9k\ 



(9) 



and find some state \(j)k) within C which has a high 
overlap with gk- 

The only aim is to find a state with high overlap, 
one does not need a certified optimal solution of 
the maximization in Eq. ([9]). Also one may replace 
in Eq. ([9|) the matrix gk by ^/gk, which may im- 
prove the convergence properties of the algorithm 
(see Section VII for more discussion) . 



3. Find an ej, > such that 



1 



9k+i ■■= 



-{gk-£k\(l)k){(l)k\) 



(10) 



has no negative eigenvalues and that furthermore 
tr(6'fc) > tr(6'fc+i) holds. In fact, one can directly 
calculate the optimal Sk, such that tr(gi|_|_j^) is min- 
imal (see Section VI); this choice is, however, not 
mandatory. In practical implementations, it can 
be useful to set an upper bound ek < emax, as this 
prevents the algorithm from subtracting too much 
from gk- 

The main idea is that if \(j)k) has a high overlap 
with gk, then it also has a high overlap with the 



eigenvector corresponding to the maximal eigen- 
value Amax(£'fc) of gk- The construction of gk+i 
leads typically to Amax(£'A:+i) < \-nax{gk) and, due 
to the normalization, Amin(£'fc+i) > Aniin(£'fc) holds. 
Hence, gk+i will be closer to the maximally mixed 
state than gk and is more likely to obey conditions 
as in Eq. ^. 

4. Check, whether gk+i fulfills conditions like Eqs. ([71 
[5|). If this is the case, then gk+i is separable and 
due to Eq. Q also gk and finally pi are in conv(C). 
Then, the algorithm can terminate. 

5. If gk+i does not fulfill Eqs. ^ [8]) return to step 2 
and k 1-^ k + 1 and iterate further until Eqs. ([71 [5]) 
hold for some k. 



Before discussing and extending this algorithm in de- 
tail, two facts must be mentioned: First, it is of course 
not guaranteed that for a given state in conv(C) the algo- 
rithm will terminate after a finite number of steps. So we 
do not claim that the algorithm can in general prove that 
a state g is in conv(C), we only claim (and demonstrate 
in this paper) that the algorithm is a powerful tool which 
works very well in practice. 

Second, a crucial step in the algorithm is the optimiza- 
tion in Eq. 1^. As already mentioned, one does not need 
a certified solution, but still it is important to find a good 
approximate solution. Clearly, the difficulty of this task 
depends on the structure of C. 

For the simple case that C are the pure bipartite prod- 
uct states, one can do this as follows: For the optimal 
10) = |a)|^) the part \a) is the eigenvector corresponding 
to the maximal eigenvalue of Xa = ti B{9k^ <8) |^)(^|) and 
\b) is similarly the vector corresponding to the maximal 
eigenvalue of Xb = tTA{gk\a}{a\ ^ !)■ This can be used 
to tackle the maximization iteratively: Starting from a 
random \a) one computes the optimal \b) via Xb, then 
with this \b) the optimal |a'), then again the optimal \b') 
etc. In practice, this converges quickly against the de- 
sired solution. For multiparticle fully separable states, 
this can be done similarly J22| . 

If C denotes the SLOCC equivalence class of some pure 
state, however, it is not so clear how to perform the 
optimization in Eq. ©. For pure three-qubit W class 
states one may use the explicit parameterization of pure 
W states from Ref. [l3l, but for more qubits, such explicit 
formulae are not available. A central step to extend the 
algorithm from Ref. Q to SLOCC classes is therefore a 
simple algorithm for the maximization in Eq. ([9|). Such 
an algorithm will be described in the next section. 

Finally, note that the termination conditions Eqs. ([71 
[S]) can also be used for SLOCC classes: From any 
pure state, one can obtain all pure product states by 
(non-invertiblc) SLOCC, hence the fully separable mixed 
states are a subset of conv(C). 



IV. MAXIMIZING THE OVERLAP OF 
n-PARTITE STATES VIA SLOCC 

As mentioned in the previous section, a crucial part of 
the algorithm is to perform the maximization in Eq. ^. 
If C is the SLOCC orbit of a suitably chosen n-partite 
pure quantum state |$o) in H = ^f^iT^i, the state after 
a general SLOCC operation is given by 



1$') 



=i^d^o 



y^($oi[®r=i4][®r=i4:]i<5'a 



(11) 



where Ai is the local filtering (or SLOCC) operator of 
the i-th party. 

The goal is then to maximize the overlap of |$') with 
a given quantum state g by applying such an SLOCC 
operation, so one has to compute 



max($'| p\^') 



(12) 



The general optimization over the tensor product of 
SLOCC operators is a hard task. Therefore, one may 
consider an iterative procedure where in each iteration 
step (i^) we optimize the overlap (fidelity) with respect 
to a single party (i) i.e., we calculate a new state by 
applying a local SLOCC operator of the i-th party and 
using the identity for the remaining parties, 



l<i>^> = 



1, 



\i 



'A^\<^^^i) 



\/<'^^-l|lnV. 
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(13) 
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here ln\i denotes the identity operator on all parties 
except the i-th party. In the following, we will usually 
omit the symbol 1„\ ^ , when there is no risk of confusion. 

In each iteration step the optimizing party is changed, 
e.g. by going from the first to the second up to the n- 
th party and then starting with the first party again. 
This iterative procedure is continued up to a fixed point, 
where the state does not change anymore. The calculated 
SLOCC operator Ai is in this case up to some factor 
proportional to the identity. Note that in general this 
optimization may converge to a local cxtremum only, but, 
as discussed above, global optimality is not required for 
the separability algorithm. 

A possible way to deal with Eq. ([T3|) is to perform a 
direct numerical optimization over the Ai . If di is the 
dimension of the local Hilbert space Hi this requires an 
optimization over 2d^ — 1 real parameters. For multi- 
qubit states this is directly feasible, however, for larger 
local dimensions or for a large number of particles n, it 
is necessary to have an analytical method to find an Ai 
which increases the overlap. This analytical approach 
will be explained in the following. 

Increasing the overlap in each iteration step is equiva- 
lent to (for brevity we use i = iip) 



(*vle|$. 
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(14) 



or equivalently 

- ($^_i| g |$^_i) ($^_i| 4A |$^_i) > 0. (15) 

In this inequality only the operator Ai is unknown. 
We denote the overlap (fidelity) of the previous iteration 
step by F^_i = ($y_i| p |<i><p_i) and choosing a local 
orthonormal basis we can rewrite 



A = ^ar,s\r){s\, 



(16) 



9 = X! '''■«i.sC2kj6)(s,6l and (17) 
\%-i) = Y.''ri\r,0, (18) 

with Qr^s e €, rr^i,si2 = Ki2,ril ^ ^' ^^'^ '^''« ^ ^' 

where we used the multi-indices ^i , ^2 which denote all 
index elements of the n \ i-partitc system. 

With this parametrization the last term of the left hand 
side in Eq. (|T5|) takes the form 



(*^-i|AIA|$^-i) = tr(4Atr„v(l<I'^-i>(*^-il) 
- tr(4Ac) =Y.a*i^ai^{C)hj 



h.j,i 

'IJV Cl,C2 

= {a\C\a). (19) 

Here, C is just the matrix representation of 
tr„\i(|$,^-i)($,p-i|) and we used {C)hh,i23 = KhC'hj, 
that is 

(<5)ci=(ih),C2=(y) - (tr„v(l*^-i)(*v-il)),y 

{C)^-^=(^ih.)X2=(^^ij) = (20) 

and (^1 , C2 is a mapping of a two valued index (the "ma- 
trix element" indices of the i-th party) to a single valued 
integer index. 

The remaining unknown expectation value of the left 
hand side in Eq. (fT5|) becomes 



(^^.ilAjpA, |$^_i) = Y^ a;zqV^''Cij?2"j™^™«2 
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(21) 



where 



(-D)ci=(ftOX2=(jm) = z2 ^ki'^Ht,j^2'^m^2- (22) 

?i,e2 



Therefore, we can rewrite Eq. (jlSp as 
(a\D^F^^xC\a) > 0, 



(23) 



where F^_i = ($^_i| £< |$„,-i 



Note that the matri- 



ces D and C are hermitian and F^p^i is non- negative. 
The maximum left-hand side corresponds to the maxi- 
mal eigenvalue (Amax) of the matrix D — F^-iC or like- 
wise the left hand side of Eq. (j23p is maximized by using 
the eigenvector corresponding to the maximal eigenvalue 
(Iflmax))- By undoing the mapping we obtain the SLOCC 
operator, 



A. 



(24) 



This procedure gives the following insight into the opti- 
mization over one local filter as in Eq. (jl3[) : If the max- 
imal eigenvalue of D — F^p-iC is positive, one can still 
increase the overlap with a suitable Ai. The correspond- 
ing eigenvector gives an Ai which increases the overlap. 
Note, however, that this Ai is optimal for Eq. (fT5|) . but 
not necessarily the optimal A^ for (4>^| g |$^) in Ec^. [TH 
In the practical implementation, especially at the be- 
ginning of the optimization procedure it is not helpful to 
use Ai directly as the SLOCC operator, because this op- 
erator is not necessarily invertible. This could therefore 
correspond to an irreversible operation which destroys 
entanglement. It turns out that using 



J±i — J±i -\- Ai] 



cl,: 



(25) 



avoids this problem. With this SLOCC operator the 
state for the next iteration step is calculated according to 
Eq. (J13p . During the iteration procedure, Amax will de- 
crease [see Eq. (fT5|) ]. the SLOCC operators become close 
to the identity, and the convergence criterion is that Amax 
is up to numerical precision zero. 



V. FINDING THE OPTIMAL e^ 

A second optimization occurring in the algorithm is the 
task to find the best Sk (see Eq. [T0|) . In detail, we want 
to maximize the decrease in the purity in each iteration, 
that is 



max 



'_tT{gl) - ti- {gl+,)] 



(26) 



with gk+i = (e^fc -efc|0fc)(</'fc|)/(l -Sfc)- With the ab- 
breviation {(pkl Qk \<t'k) — c the above maximization leads 
to 



max 



4[tr(gD-l]+2gfc[c-tr(e2,)] 



(27) 



Taking the derivative with respect to £/c, we find the max- 
imum as 



c- tr(gg) 



1-c 



(28) 



In the implementation, in case of e™'^^ > 10 ^A^, we 



define e? 



10 Xd, where A^ is the minimal eigen- 



value of gk (guided by practical experience), to keep the 
remaining state positive during the iterations. This cor- 
responds to an upper bound on e^, as mentioned above. 



VI. INCREASING THE MIXEDNESS 

A central strategy of the algorithm is to increase the 
"mixedness" of gk in each iteration step i.e., to lower 
the purity tr(g^) . One may wonder, whether this is al- 
ways possible by subtracting some \4>k){(l>k\ S C. Indeed, 
one can show that this is the case, unless gk is outside 
conv(C), which means that the algorithm has no chance 
to succeed anyway. 

To see this, the condition for an increase of the mixed- 
ness can be formulated as 



ti- (gl) > tT{gl+,] 
1 



tr 



Qk - £k\(t)k){(t>k\ 



{i-Sk)' 
1 



I- £k 
[tr {gl) - 2ektT{gk\cj)k){c^k\) + el] 

tr{gl)-2ek{(l)k\gk\<Pk)] 



l-2ek 
for small Sk- It follows that iff 

{<l>k\gk |0fc) > tr(gfe) 



(29) 



the state \(j)k){4'k\ can be subtracted with a small weight, 
and the mixedness increases. 

If condition Eq. ((29l) is not fulfilled for any |0fc) the 
mixedness cannot increase. But this implies that 



sup ' 

I0>ec 



\Qk 



< tr(g^) 



Consequently, tr(g^.yV) < for the observable 
W^al-gk 



(30) 



(31) 



with a = sup|j,\gg (01 gk \4>) ■ This W is nothing but a 
witness [3| which discriminates between conv(C) and the 
remaining states, and tr(pyV) < implies that a state 
is not in conv(C). Therefore, states which cannot fulfill 
the condition in Eq. ()29p are either not in conv(C) or (in 
case that tr{gkW) = 0) they may lie at the border of 
conv(C). This, however, is a set of measure zero and not 
of practical relevance. 



VII. CONVERGENCE BEHAVIOR OF THE 
ALGORITHM 

Let us now discuss some practical issues. The question, 
whether or not the algorithm converges depends first on 
the type of state and decomposition to be determined and 
second on the distance of the state from the boundary of 



the considered convex set. The closer the state is to a 
boundary the slower is the convergence. The algorithm 
does e.g. not work with rank deficit states, because over- 
lap of the optimized pure states with the kernel of the 
density operator cannot be avoided i.e., it is not possible 
to ensure pk+i > 0. In the three- and four-qubit case the 
decompositions consists of the order of 10'^ states (mean- 
ing that the algorithm requires this number of iterations, 
until the conditions in Eqs. ^ [S]) apply) i.e., usually 
such decompositions contain many more states then the 
Caratheodory-bound of cf . 

In practice, the overlap optimization in Eq. ^ with 
the square root y/g instead of g has a better convergence 
behavior. Also other fractional powers of g show a similar 
advantage. Note that replacing g by ^^g does not affect 
the proof that the iterated state is separable, if Eqs. (0 
E) apply. 



VIII. LOWER BOUND ON THE e-BALL VIA 
CROSS POLYTOPE 

The presented algorithmic procedure allows to deter- 
mine a decomposition of the state g with the specified 
SLOCC properties. After n iterations we have a decom- 
position i.e., the set S ~ {{\4'i) {4'i\} ^ 9n} , of our initial 
state of the form 



g = '^Pi\(l)i){<l>i 



qnQn 



where the probabilities pi are given by pi 

i 

1i = n (1 -Pj) with qa := 1 and go := g. 



£iqi 



(32) 



and 



By construction, our convex set S has specific "entan- 
glement" properties which are valid for all states in its 
convex hull. However, if we obtained the state g from 
experimental data, we have to deal with errors and im- 
perfections, and so the starting state g is affected by un- 
certainties. Therefore it is of great importance to give 
some statements about the "stability" of the determined 
decomposition, or an estimate of the probability that an 
experimental states lies inside this convex set. 

A first possibility to deal with this problem was used 
in Ref. Q: There, starting from the experimentally ob- 
tained state gexp the measurements were simulated via 
a Monte-Carlo simulation, and 200 sampled states were 
reconstructed via a maximum likelihood approximation. 
Then, separability of the state gexp was only claimed, 
when the algorithm could prove that gexp as well as all 
samples were separable. Note that the generation of 
states via Monte-Carlo simulation of the measurements 
is a standard technique to estimate errors in ion-trap ex- 
periments. 

A different possibility can be obtained by answering 
the question, how much can an experimental state devi- 
ate from g such that the state still belongs to the set S. 
Although this question is in general not easy to answer 
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FIG. 1: (Color online) Schematic two-dimensional example 
for a decomposition of g (red dot). The convex hull of the 
states of the decomposition i.e., the set S, is denoted by the 
dashed line. The cross polytope is shown in blue, the vertex 
states of the cross polytope are marked by blue triangles and 
the Ebaii corresponds to the circle. 



[23|, we can determine a lower bound on the minimal 
Hilbert-Schmidt-distance of the state g with respect to 
the convex hull of S. 

The idea is to show that if the state deviates from 
g in different directions, then it remains in the convex 
set. More precisely, we construct a symmetric cross poly- 
tope [IJI with the state g in the center (see Fig. [T] for 
a 2-dimensional example). For a quantum state g with 
Hilbert space dimension d the set needs at least (P con- 
stituents s.t. a nonzero volume object in this space is 
possible. The cross polytope is a symmetric polytope 
with 2{(P — 1) vertex states. The vertex states of the 
cross polytope are defined by (f^^ = gzL fcpCi, where we 
used a vector representation of the density operator in 
((P — l)-dimensional Euclidian space by mapping 



/ai ad + iad+i 



A 



asd-s + ia3d-2\ 



ad-i 



d-l 



ad2_2 - iarf2_i 1 - J2i=i Oi / 



A: 



a2 



ya-d^-i/ 



The elements of the basisvectors in this notation are given 
by (e.)j = Sij. 

The Euclidian distance of the vertex states of the cross 
polytope with respect to the state g, which is given by 
/cp, is maximized under the constraint that the vertex 
states arc contained in the convex hull of our set S (see 
Fig. [J). 

First, using a divide and conquer algorithm we calcu- 
late the maximal parameter / of each state g±i = g± fci 
such that it is contained in the convex hull of our set S. 
Whether a state is contained in the convex set can be 
decided by using a linear program, e.g. via the Matlab 



routine linprog. The vertex state with the smallest pa- 
rameter /cp is used for defining the vertex states of the 
cross polytope inside our convex hull. The parameter /cp 
also depends on the relative orientation of the chosen or- 
thogonal basis {e,;}, here an additional optimization is 
possible. 

Then, the smallest Hilbert-Schmidt distance of g with 
respect to any point in the convex hull of the cross poly- 
tope is given by £baii 



/c, 



/-p — -, where d is the dimension 
of the Hilbert space and /cp is the maximal parameter 
such that the cross polytope is contained inside the set 
S [lj|. We will present an example below. 



IX. EXAMPLES 

In this section, we present several examples for the 
application of our algorithm outlined above. 
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FIG. 2: (Color online) The lower bound on the size of the 
£-baIl of the state given in Eq. (|33|) as a function of p. The 
decompositions are obtained with the procedure of Sec. lVIIIl 
The line is a polynomial fit of the points for guiding the eye. 
The ball-size strongly depends on the number of constituents 
in the decomposition, especially if states far away from the 
border are considered. Therefore here decompositions are 
compared which all consists of about 4000 pure states. 



A. GHZ states affected by white noise 

First, we consider the GHZ state of n qubits affected 
by white noise 



eGHZ„b)=p|GHZ„)(GHZ„ 



1-p 



1, 



(33) 



where |GHZ„ 



i5(|0...0) + 11...1)). For three 



V2 



qubits, these states have the following properties: 

• PGHZ3 is fully separable ifFp < 1/5 = 0.2 |24| . 

• PGHZ3 is genuine multipartite (tripartite) entangled 
iff p > 3/7 « 0.4286 [H, 

• eoHZs belongs to the GHZ class iS p> 0.6955. For 
a detailed discussion see Ref. [l^ ■ 

With our algorithm we can determine separable and 
biseparable decompositions up to the threshold values of 
p. A bound with a threshold value of p = ||| « 0.6944 
for existing W-class decomposition was obtained with our 
algorithm, which is close to the optimal threshold value 
ofpw 0.6955. 



1. The £-ball and the robustness depending on p 

With our algorithm it is possible to obtain W- 
decompositions of the state pghzs up to p ~ 0.6944. Due 
to the "small" number of vertex states forming the con- 
vex set the size of an e-ball generated via the procedure 
of Sec. IVIIII will be small in comparison to the maxi- 
mum possible ball which fulfills the considered proper- 
ties. However the lower bound on the size of this ball 
will also depend on the distance of the considered state 
to the border where the properties are not fulfilled any- 
more. In Fig. [2] the lower bound on the size of the e-ball 



is plotted versus the parameter p. As expected for the 
almost maximally mixed state (small p) the ball is quite 
large, but it decreases by several orders of magnitude as p 
approaches the threshold value of the considered convex 
set. 

From our e-ball we can also deduce that any state g' 
which has an Euclidean distance from g which is smaller 
equal the radius of the e-ball is contained in our convex 
set i.e., using the Hilbert-Schmidt-distance between two 
operators A and B given by d = y/ ti{{A — B)'^), we can 
calculate the distance between two GHZ- Werner states, 
see Eq. ((33|) . with purity p and p' = p + S to be given by 



d' = tr([eGHzb)-gGHz(p')]') 



(34) 



If we have an e-ball in the surrounding of gcuz (p) with 
radius Ebaii then we can deduce that also the state 
QGHzip') with 



(35) 



S < -\/-ebaii 



belongs to the same convex set, allowing to increase the 
threshold parameter accordingly. 



2. Entanglement properties o/pghz4(p) 

The maximal parameter popt ~ 0.467 for biscparability 
of the four-qubit GHZ-state mixed with white noise was 
derived in Ref. [25|. Our algorithm is able to determine 
decompositions for values of p up to p « 0.466. 

For four qubits there is already a continuous set of 
inequivalent SLOGG entanglement classes, and it can 
happen that a state can be decomposed into bisepara- 
ble states, but not into SLOGG equivalents of some gen- 
uine multipartite entangled states. To investigate this. 



we considered the state pcnziip) and asked when it can 
be decomposed into SLOCC equivalents of the four-qubit 
W-state, IWi) = i(|0001) + |0010) + |0100) + |1000)). This 
seems to be only possible for p < 0.32, but for the larger 
parameter regime p < 0.467 the state is biseparable. 

An intuitive argument for such a behavior is that 
the four qubit W-state can be transformed via SLOCC 
into a thrce-qubit GHZ state as well as to a two-qubit 
Bell state tensored with a product state, but it is not 
possible to reach |04) = |0+) (8)1^+), we even have 



i-na^WGSLOCCiWi) \{W\(I)4)\'^ = 0.5[3J|. States like |04), 
however, are essential in the biseparable decomposition 
of PGHZ4 (p) 123 . Many lower entangled pure states are 
SLOCC inequivalent to specific genuine multipartite en- 
tangled pure states. 



B. W-states with white noise 

In order to give an example where the algorithm is not 
capable of computing the threshold of separability, we 
consider states of the form 



QWnip) =p\Wn){Wn 



l-P 

2" 



1, 



(36) 



where the three-qubits W state \W^) is given in Eq. ^ 
and for four qubits we have IW4) = ^(lOOOl) + |0010) + 
1 0100) + 1 1000)). The border to the class of bisepara- 
ble states can only be roughly approached by the algo- 
rithm. In the three-qubit case we have a gap of about 
Ap = 0.03, the exact value p = 0.4790 is known from 
Ref. |2a]. For the four-qubit W-state with white noise 
the gap with respect to the upper bound obtained by 
an semidefinit-programming-witness (SDP-witness) [2Q| 
is approximately Ap4 = 0.04 i.e., we determined a de- 
composition for p = 1 - 0.526 - 0.04 = 0.434. 



C. Bound entangled state from an unextendible 
product basis 

As a second example, we consider the bound entangled 
states arising from an unextendible product basis [27|. 
These states are defined via using the product vectors 

1^0) = |0)(|0) - |1))/V2, iV'i) = (|0) - |1))|2)/A IV'2) = 
|2)(|1) - |2))/V2, IV3) - (|1) - |2))|0)/V2, 1^4) = (|0) + 
|1) + |2))(|0) + |1) + |2))/3. Then the state 



QBE, 



=\U-ii\^'^){w 



(37) 



is an entangled state on a 3 x 3 system, which is not 
detected by the PPT criterion. We considered the family 
of states 



£'upb(p) = PQbf. + (1 - p)l/9. 



(38) 



They have often been used as a test-bed for separability 
criteria. To our knowledge, the best criterion for these 



states is the first step of the algorithm of Doherty et al. 
j28| which detects them to be entangled for p > 0.8691. 
Our algorithm proves that these states are separable for 
p < 0.83. 



D. Thermal states wth the Heisenberg interaction 



Let us consider the thermal state 



gniT) - cxp{-Hh/T} 



(39) 



of three spin-1/2 particles interacting with the Heisen- 
berg interaction. 



Hh = / hij with hij = 



k—x,y,z 



.(*) 



M) 



(40) 



where i,j e {1,2,3}. In Ref. [231 the entanglement prop- 
erties of this system were investigated, and it was shown 
that the spin-squeezing inequality (AJ^.)^ + (AJy)^ + 

(AJ;,)2 > N/2 with Jk = 1/2 Ei <^k^ detects these states 
as entangled for T < 5.461. Remarkably, the spin- 
squeezing inequality shows that for 4.329 < T < 5.461 
the thermal state is biseparable with respect to any bi- 
partition, but not fully separable. Direct application of 
our algorithm gives that for T > 5.462 the thermal state 
is fully separable, giving strong evidence that the spin- 
squeezing inequality is a necessary and sufficient criterion 
for the thermal state. For more than three spins, how- 
ever, this does not seem to be the case. 



E. Experimental pseudo bound entangled state 

In Ref. [ij] a class of three-qubit bound entangled 
states with rank seven were introduced. We consider 
the specific state where the entanglement is maximally 
robust with respect to white noise |30[ : 



PBE3 



1 

TV 



2|GHZ)(GHZ|- 



a|001)(001| -t- a|010)(010| + -|011)(011k 

a 

a|100)(100| + -|101)(101| + -|110)(110| 
a a 



(41) 



where |GHZ) = -^ (|000) + |111)), a = 0.3460, and the 

normalization is A^ = (2 -f- 3(a -I- i)). The state has the 
curious property that it is biseparable with respect to 
any bipartite splitting, but it is nevertheless entangled. 
A pseudo bound entangled state of this form was experi- 
mentally generated and characterized in Ref. [3i|, see Fig. 
[21 It was shown via a witness operator that the state is 
entangled and PPT with respect to any bipartite split- 
ting. Now with our algorithm we are able to prove that 
this state is biseparable with respect to the split B-AC 
with an e-ball of ebaii = 4 • 10""'. 
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FIG. 3: (Color online) Experimental pseudo bound entangled 
state [31( . All imaginary elements of the experimental state 
are small and therefore not shown. 




column index 



FIG. 4: (Color online) Real part of the experimental W-state 
from Roos et al. 132] (left) in comparison to the ideal W-state 
(right). All imaginary elements of the experimental state are 
small and therefore not shown. 



G. Summary of the algorithmic performance 

In this section we summarize the obtained threshold 
parameters for the various examples considered. Our 
standard approach was to consider different types of 
quantum states mixed with white noise of the form: 



Qi{v) =VQi + (1 -]5)l/d. 



(42) 



Then we determined the threshold parameters p for var- 
ious types of entanglement which we can achieve via our 
algorithm, and to compared them with bounds from the 
literature or with exact values if they are known. The 
only exceptional parameterization was the case of the 
thermal equilibrium Heiscnberg spin chain state (Sec. 
lIXPp which depends on the temperature. In Table U 
we summarize the threshold parameters reached via our 
algorithm and show the best known bounds from the lit- 
erature. For several cases we also calculated lower bounds 
on the e-ball for our determined decompositions (see Sec. 
IVIIip . For all cases where the exact bounds are known 
we can reproduce the threshold parameters quite well 
by our procedure, apart from the three-qubit W state, 
where there is a small gap. It is remarkable that our 
algorithm seems to work independently of the type of 
decomposition to be determined. For the given exam- 
ples the calculation needs less than 5 min on a standard 
personal computer. The number of terms calculated in 
the decomposition strongly depends on the distance to 
the border of the convex set. I.e. already a very small 
change on the p-threshold values has a huge influence on 
the number of terms in the decomposition. For threshold 
values with the given precision in Table U it is possible to 
find decomposition within the order of ~ 5000 terms. 



X. CONCLUSIONS 



F. Experimental three-qubit W-state 



Nowadays it is possible in several experimental se- 
tups to generate quantum states which are e.g. close 
to the three-qubit W-state. The generated states are 
characterized by quantum state tomography. Here we 
use data from an experiment by Roos et al. [32|. This 
"typical" experiment (sec Fig. |4]) led to a fidelity of 
F = \/{W\ Qcxp \W) =0.9, where we can prove via an 
entanglement witness that this state is genuine multi- 
partite entangled. Now the question arises whether this 
state really belongs to the W-class of entanglement. With 
our algorithm it is not possible to find a W-class decom- 
position of goxp, so probably the state belongs to the 
class of GHZ-entangled-states. On the other hand we 
can find for the slightly depolarized version of this state, 
e.g. g' — 0.96pcxp + ^^^l, a W-dccomposition and we can 
prove that it is genuine multipartite entangled. 



We presented an easy to implement, fast and straight- 
forward method for finding decompositions of quantum 
states with specific SLOCC entanglement properties. For 
a large variety of examples decompositions were deter- 
mined for those parameters, where they are known to 
exist. For instance, is it possible to find separable de- 
compositions of Werner states with minimal amount of 
white noise. Also the threshold values for biseparability 
can be reproduced in most cases for three- and four-qubit 
states. Especially interesting are the cases where the ex- 
act parameter range for the existing decompositions is 
not known. It was e.g. possible to investigate bound en- 
tangled states with our algorithm. This method of deter- 
mining specific decompositions of a quantum state is a 
complementary tool to entanglement criteria such as en- 
tanglement witnesses. Together these tools allow to ex- 
tensively specify the entanglement properties of a given 
quantum state. In the future we would like to understand 
better the convergence behavior of the algorithm, espe- 
cially why optimally decreasing the purity of the quan- 
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State 


Ent. 


Bound 


Decomp. 


£-ball 


PGHZ3 


S 

BS 
W 


1/5* (a) 
0.429* [26] 
0.6955* fl3] 


0.199 

0.4285 

0.694 


8.1 • 10"^ 

9.2 • 10"'' 
7.1 • lO-'^ 


eGHZ4 


s 

BS 

W 


1/9* (a) 
0.467* [26] 


0.111 
0.466 
0.316 




ew3 


s 

BS 


3/11 (a) 
0.479* [26] 


0.1727 
0.45 


7.6 • 10-* 
1.1- 10"^ 


PW4 


S 
BS 


1/5 (a) 
0.474 [26] 


0.09 
0.434 




PUPB 


s 


0.87 [28] 


0.83 


1.2 • 10"'' 


QBES 


S 
BS (AB-C) 


0.786* [33] 
1* [12] 


0.726 
0.9 


2.0 • 10"'' 
1.1- 10"^ 


QHiT) 


S 
BS (AB-C) 


5.61 [2^ 
4.329 [29] 


5.62 
4.33 


2.1 • 10"^ 

7.2 • lO"** 



turn state in each iteration is at the end not necessarily 
a good strategy. 



TABLE I: Threshold values for p, see Eq. (gg]). With an aster- 
isk we denote exact values from the literature and with (a) we 
denote bounds obtained via the PPT criterion. The column 
Decomp. contains the parameter threshold up to which we 
are able to determine S (separable), BS (biseparable), W (W- 
type) decompositions. The specific methods used to obtain 
bounds on the threshold values are discussed in the corre- 
sponding paragraph of the examples. 
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